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Abstract 

When one slightly pushes a thin elastic sheet at its center into a hollow cylinder, the sheet 
forms (to a high degree of approximation) a developable cone, or d-cone for short. Here we 
investigate one particular aspect of d-cones, namely the scaling of the elastic energy with the 
sheet thickness h. Following recent work of Brandman, Kohn and Nguyen [2] we study the 
Dirichlet problem of finding the configuration of minimal elastic energy when the boundary 
values are given by an exact d-cone. We improve their result for the energy scaling. In 
particular, we show that the deviation from the logarithmic energy scaling is bounded by a 
constant times the double logarithm of the thickness. 

1 Introduction 

Stress and energy focusing in thin elastic sheets has recently attracted a lot of interest in the 
physics literature [10, 3, 5, 12, 6, 8, 9, 11, 4]. One basic feature are (almost) conical singularities. 
A conical singularity arises, e.g., in the following experiment. Put an elastic sheet of radius 1 
concentrically on top of a hollow cylinder of radius R < 1 and push the sheet down at its centre. 
It has been observed that the sheet assumes (to a high degree of approximation) the shape of a 
developable cone (or d-cone for short). In the physics literature, this has been discussed e.g. in 
[3, 5, 10, 12]. There are several remarkable features of the d-cone: The angle subtended by the 
region where the sheet lifts off the rim of the container is a universal constant (approx. 139°), 
independent of the indentation, the thickness and the material of the sheet (for small indentations, 
[5]). The tip of the d-cone consists of a crescent-shaped ridge where curvature and elastic stress 
focus. In numerical simulations it was found that the radius of the crescent i? C res scales with the 
thickness of the sheet h and the radius of the container i? C ont. as R C i CS . ~ ^^Rcint ■ This depen- 
dence on the container radius of the shape of the region near the tip is not understood [12]. As 
argued in this latter reference, it cannot be explained by an analysis of the dominant contributions 
to the elastic energy, which are: The bending energy from the region far away from the center, 
that is well captured by modeling the d-cone as a developable surface there; and the bending and 
stretching energy part from a core region of size 0(h) where elastic strain is not negligible. The 
result of this (non-rigorous) argument is an energy scaling E ~ h 2 (Ci \ \ogh\ + C2). 
Here, we discuss the scaling of the elastic energy with ft in a rigorous setting. 



The natural variational formulation is to minimize the elastic energy of the sheet (which contains 
stretching and bending contributions) in the class of deformations y : B\ — > R 3 which satisfy the 
obstacle constraint 

Im(y) n {x G R 3 : x\ + x\ = R 2 cont ,,x 3 < H cont .,y(0) = 0} = 
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where Bi = {x E R 2 : \x\ < 1} is the reference configuration of the 2-dimensional sheet, 
and -ffcont. is the height of the container. The derivation of a precise lower bound of the form 
E ~ h 2 (Ci \ \ogh\ + C 2 ) for this model looks very hard. Indeed even, for the much more severe 
constraint that y maps B\ into a very small ball the only known rigorous lower bound is that 
lim/ l _ ! . 00 Eh(y)/h 2 = oo while it is conjectured that the correct scaling is Eh(y) ~ Ch 5 ^ 3 in this 
setting (see [7] for a discussion and a rigorous proof of the upper bound). To make progress in 
rigorously understanding the asymptotic influence of the regularizing effect of the thickness h we 
follow recent work of Brandman, Kohn and Nguyen [2] and free ourselves from the specific obstacle 
type constraint and consider instead general Dirichlet problems where the boundary conditions are 
given by an exact developable cone. 

This also partly motivated by results that have been obtained in the physics literature for 
a simpler model. In [5], the class of allowed deformations y is restricted to isometrics with a 
singularity at the origin. These maps are completely determined by their values on the boundary 
dB\. In order for these maps to be isometric away from the origin, the boundary values have to 
be unit speed curves 

y\ dB = 1 :dB 1 ^S 2 . 

This effectively reduces the problem from a 2-dimensional to a 1-dimensional one. After a suitable 
ad- hoc renormalization of the bending energy (i.e., cutting out a small ball around the origin where 
the energy density becomes singular), the elastic energy is minimized as a functional of 7. In this 
simpler setting, an obstacle of the above type is treatable. In fact, the shape of 7 in the region 
where the sheet lifts off the rim of the cylinder is prescribed by an ODE that can be solved more 
or less explicitly (for small indentations, see [5]). 

2 Setting and statement of the main theorem 



Let B r = {x e R 2 : \x\ < r}, A r = B r \ B r/2 and for y e W 2 ' 2 {B U R 3 ) let 

Eh(y)= I (\Vy T Vy-Id\ 2 + h 2 \V 2 y\ 2 )dx. (1) 

J Bi 

Furthermore, for a curve 7 <E C 3 (dBi, R 3 ) with = | = 1, let 

V 1 = {y G W 2 - 2 {B U R 3 ) : y\ dBl = 7, y(0) = 0} . 

In the following, consider such a 7 to be fixed. By y(x) — \x\j(x), we denote the 1-homogeneous 
surface with boundary values prescribed by 7. 

Existence of minimizers of Ef t in the class V 7 follows easily from the fact that Eh is coercive 
and convex in the highest derivatives and the compact embedding W 2 ' 2 c — > W 1 ' 4 . Our main result 
is 

Theorem 1. Suppose that 7 does not lie in a plane. Then for sufficiently small h we have 

Ci In 1 - d In fin ^) - C 2 < ^ min E h (y) < d In \ + C 3 , 
h \ h J h z yev-, h 

where 

Ci=CT 1 ( 7 ) = -i- f \V 2 y\ 2 dx 
In 2 J Bl \B 1/2 

and C2, C3 only depend on 7. 



This improves work of Brandman, Kohn and Nguyen [2] who showed that 
liminf fe ->o /^m(i/fe) minE h > Ci/2 and limsup, t ^ h i ln \ 1/h) minE h < C\ and, very recently in 
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parallel to our work, that lim^o h2 ^77^ rnin Eh — C\ . The result above shows that the deviation 
from the leading order logarithm is at most a double logarithm. Indeed a natural conjecture is 
that this error is or order 1, but we have not been able to prove or disprove this so far. 

Our proof of the crucial lower bound consists of three steps: First we estimate the L°°-norm 
of y in the ball Bh (sec Lemma 4 below). Then we derive the key estimate for the L 2 -norm of 
e = y — y on dyadic rings A 2 -j (see Lemma 5 below). Finally we argue that since y is close to y in 
L 2 , the bending energy of y can be bounded from below by the bending energy of y, up to a small 
error. 



3 Upper bound and L 2 estimate of e = y — y 

The proof of the upper bound is standard and we include it for the convenience of the reader. 
Lemma 1. 

inf E h (y) <C 1 h 2 \n\ + C 2 h 2 

yev~, h 



W / le r e C 1 = (ln2)- 1 / BABi/2 |V 2 y| 2 . 

Proof. Let : R -> R be a C 2 function with <j>(t) = t for t > 1 and (j){t) = for t < 1/2 and set 

y h (x) = H{j^) 7 
Then yh = y on B\ \ B] x . Hence (Vy) T S7y = Id in B\ \ Bh and 

f (|Vy T Vy-/rf| 2 + ^ 2 |V 2 2 ;| 2 )dx-C 1 lni (2) 

JBi\-B h h 

Moreover in Bh one has the estimates \Vyh\ < C and |V 2 y/j| < C/h. This implies the assertion. 

□ 

Now we prove some auxiliary lemmas which will allow us to estimate sup Bh \y\. This estimate 
will be needed in the proof of the L 2 bound for of y — y on dyadic rings. 

Lemma 2. Let v e W 2 - 2 {B h ). Then 



sup 

x£B h 



v(x) - t>(0) - ( / Vv(x')dx' 

\JB h (0) 



< Ch\\V 2 v\\ LHBh) 



Proof. For h = 1 this follows from the embedding W 2 ' 2 (Bi) ^ C°(Bi) and the Poincare inequality. 
For h 7^ 1 the assertion follows by considering the rescaled function u{x) — j^v(hx). 

Lemma 3. Let w e W 1 ' 2 {B 1 ), < e < 1. Then 



□ 



-r wdx — -f 

\JB e JB X 



w dx 



< C ( In - 



1/2 



IVwI dx 



1/2 



Proof. The Poincare inequality for w — w implies that for R g [1/2, 1) 



-r wdx — -f wdx < C 

JB R J B 1 JE 



< C I \Vw\dx 

Bt 



(3) 



and scaling yields in particular 



-/ wdx—-f wdx < C- [ \Vw\dx 

J B r/2 J B r r J B r 
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Apply this with r = 2 k for k = 0, . . . , n — 1 and define 



fix) :-E 2 V 



(4) 



k=0 



Then 



/ wdx- -f wdx <C f f\Vw\dx<\\f\\ L 2 iBl) \\Vw\\ L 2 {Bl) 



Now we have / < 2 fc+1 in B 2 - fc \ B 2 -k-i for fc < n - 1 and / < 2™ +1 in B 2 -„. This implies 
that ||/||i,2( Bl ) < Cy/n. Choose n such that 2"™ > e > 2~ (n+1 \ Then scaling of (3) and the 
Cauchy-Schwarz inequality yield 

£ wdx-1 wdx <C2 n f \Vw\dx < C\\Vw\\ L 2 {B2 _ n) 



which completes the proof since n < ln 2 -. 

Lemma 4. There exists constants such that for all < h < 1/4 and all y € V-y we have 



sup\y\ < Ch + Ch In 



B h 



1/2 



\V 2 y\\ L 2 {Bl) <Ch + C[ln 



1/2 



El /2 iv). 



□ 



(5) 



Proof. This follows from Lemma 2 and Lemma 3 (applied with w = Vv) and the following calcu- 
lation 

Vy da; 



y<gii/dS= / j®vdS 

9Bi J9Bi 



which yields 



^3 



Vy da; 



< 



2tt 



□ 



We now come to the key estimate for the difference between a low energy map y and y in the L 2 
norm on annuli. The idea is to look at fibres in radial direction in the domain, i.e. at line segments 
connecting the origin with dBi . By the upper bound on the elastic energy, the "stretching" of y 
on such a line segment is small. The boundary values of y on the line segment are fixed as well, 
and so the deviation of y from the straight line connecting the boundary points cannot be large. 

Lemma 5. Let h be small enough, 2h < r$ < 1, and assume that E h {y) < 2Ciln±. Then 

2 



i 



1 



\y - y\ dx < Cr A h In - + Cr 2 h 2 In 

lB ro \B ro/2 h \ 

where C is a constant that only depends on 7. 

Remark Note that the second term on the right hand side of the estimate is controlled by 
the first as long as tq > h In 4 • 



Proof. We consider polar coordinates and set 

r](r,6) — y(rcos9, rsin9), e(r,8) = (y — y)(r cos 9, r sin 8) 



(6) 



and we write e' = (d/dr)e and 77' = (d/dr)n. By Fubini's theorem the maps r i-> r]{r,0) and 
r e(r, 9) are weakly differentiable for a.e. 9. 
1. The key estimate is 



/ / \e , (p,9)\ 2 dpd9<Ch\n-. 
Jo Jh h 



(7) 
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To prove this note that \r/'\ 2 = \y + e'\ 2 = 1 + 2j ■ e' + \e'\ 2 which yields 

\ef = W\ 2 -l-(2 7 .e)' 

since 7 depends only on 9 but not on r. Using the orthonormal basis x/\x\, x^/lx] we get the 
pointwise estimate 

|(Vy) T Vy-/d| 2 >(M 2 -l)) 2 . 

By the boundary condition we have e(l,6) = and Lemma 4 yields \e(h,6)\ < Chlnj^ (since 
\y(x)\ — \x\). Thus an application of the Cauchy-Schwarz inequality gives 

2tt pi 

2 



\e'(p, 9)\ 2 dpd9 

< I f (|V| 2 -l)dpd^ + C/iln^ 
Jo Jh h 

<(/ / (W\ 2 -lfpdpdO) [j I -dpd9) +Ch\n 



h 
r l 



2?r r l \ !/ 2 / r 27T r l ^ \ V 2 ^ 



Jh P J h 



( 1 \ 1/2 1 

^(y) 1 / 2 (27rln-J +Chln- 



and (7) follows from the assumption on Eh(y). 
2. The Cauchy-Schwarz inequality yields 



|e(r,0)-e(M)| <r 1/2 (J* \e'(p,9)\ 2 dp 
for a.e. 9. Taking the square, integrating over 9 and using (7) we get 

f \e(r,9)-e(h,9)\ 2 d9 < OMn i . 
Jo h 



1/2 



By Lemma 4 we have 

|e(M)| 2 < Ch 2 ( In 



r 2 

h 

and the assertion follows by integrating these two inequalities from r /2 to r with respect to the 
measure r dr. 

□ 



4 Proof of Theorem 1 

We first recall two interpolation inequalities for Sobolev functions. 

Lemma 6. Let SlcR™ be a bounded Lipschitz domain, v <E W 1,2 (Q), u e W 2 - 2 (Q). Then 

IMli^an) <C|MlL2(Q)IMIvi/i. 2 (f}), (8) 
||V«|| 2 2(n) <C||u||£ a(n) +C|M| L2(0) ||V 2 u || i2(0) (9) 

w/iere oZZ constants C only depend on SI. 

Proof. The first inequality follows from the continuity of the trace operator W 1 / 2,2 ^) — > L 2 (<9S1), 
and the fact that W /1 / 2 - 2 (Q) is a real interpolation space of the pair (L 2 (f2), W 1,2 (SV)) which yields 

\Hh{en) < C\MwW(n) ^ c \\ v h'(a)\Mw^'(a) , 
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see [1]. The second interpolation inequality follows directly from Theorem 5.2. in [1] which states 
that for all e < 1 

l|V«|| ia( n) < K (^\\u\\ L 2 {a) + e\\W 2 u\\ L 2 (n) 
If || V 2 m||l2( ) < ||u||l2(q) onecantakee = 1 to obtain (9), otherwise one takes e 2 = ||u||i,2(fi)/|| V 2 m||i,2(q) 



□ 



Proof of Theorem 1. Let M e N (to be chosen later). Recall that e = y — y. We have 
^E h (y) > [ |V 2 y| 2 dx 



J Bi 



\B 2 _ M 



|V 2 y| 2 dx 



- 2 / 

JBi 



V 2 e:V 2 ydx + f \V 2 e{ 2 dx 

\B 2 -m JBx\B 2 ^ M 



Consider the second term on the right hand side with the integration restricted to the annulus A ro . 
We define e : A\ — > R 3 by e(x) = 1 e(r x). Observe that 

V 2 e(.x) = r (V 2 e){r x), V 2 y(a;) = r (V 2 y)(r x). 

Since y is 1-homogeneous we get 

/ V 2 e:V 2 ydx= J rl(W 2 e)(r x) : (W 2 y)(r x) dx = V 2 e:V 2 ydx. 

J Ai J Ax J A TQ 

We integrate by parts to obtain 

/ Vf.Vyto-f tM . fM *-{ e,,. S , w dx + / 4 ,.*^dS, 

,/Ai JAj JAi JaAi 

where v is the unit outer normal on dA\. Hence 



/ 

J Ax 



V 2 e : V 2 yda; 



<l|V 3 y||L2( Al) ||Ve|| i 2 (yll) + ||Ve|| i 2 (aAl 

<c(||Ve|U2 (Al) + ||Ve||^ 2 2 (Ai) ||Ve||^, 2(Ai) ) 
<c(||Ve|| L 2 (Al) + ||Ve||^ 2 (Ai) ||V 2 e||^ (Ai) ) 

where we used Lemma 6. Applying again the second inequality in that lemma, we get 



/ 

J Ax 



V 2 e : V 2 2/dx 



<C lie 



(» 



L 2 (Ax) 



1/2 

\L*{Ax) 



||V 2 e 



2aiiV2 



!i 2 (A!) 



1/4 

\L?(Ax) 



||V 2 e 



2-||3/4 



(10) 



We apply Young's inequality ab < -S p a p + \8 p 'b p ' with the pairs p = |, p' = 4 and p = I, 



p' = |, and <5 p = ^ to obtain 



V 2 e : V 2 ydx 



< C (\\e\\ LHAl) + \\e\\%\ Ai) + \\e\\%\ Ai) )) + 1||V 



2 e||i2 (Al) 



(11) 



Now we undo the rescaling. We have 

\\e\\L*(Ax) = ro 2 \\e\\ L 2 (Aro) , || V 2 e|| L 2 (Al) = || V 2 e\\ L 2 {Ara) 

By Lemma 5 

^He|U^ )<c(A) 1/2 (ln^ 



1/2 



(12) 



as long as 



Thus 



ro > ft In — . 
ft 



J 

JA, 



V 2 e : V 2 ydx 



<C[- 

ro 



1/5 



-I 



1/5 



2" v2e lli^o) 



(13) 



(14) 



as long as (13) holds. 
Choose M such that 



lo S2 ^ - lo §2 In i < M < log 2 i - log 2 In i + 1 . 



(15) 



This implies that 



2 -(M-i) >Mn 



1 



In particular for sufficiently small h and for ro > 2 ( M x ) the inequality (13) holds and we also 
have ro > 2ft. We thus get 



f |V 2 y| 2 dx 

7 Bi\B 2 -m 



>[ \V 2 y\ 2 dx-2 [ V 2 y:V 2 edx+ [ |V 2 e| 2 da; 

Jb 1 \B 2 . m Jb 1 \B 2 ^m JB!\B 2 -m 

>M (j \V 2 y\ 2 dxj - 2C J2 2 ' /5hl/5 ( ln 
>M^j \V 2 y\ 2 dx^j -20C2( M ~ 1 V 5 h 1 / 5 (in^ 5 
^In2( 1 ^- ln ( 1 ^))(/. 1 |V2 ^ dX )- 20C 



This completes the proof of Theorem 1. 



(16) 



□ 



Remark. The estimate (16) in connection with the upper bound on Eh shows that for any y 
with Eh(y) < C\ In ^ + C3 (and in particular for any minimizer of Eh) we have 



j |Vy T Vy- Id\ 2 dx< Ciln^ln^ + C 4 , 
j \V 2 y\ 2 dx<G\\n(\n^j + C 4 , 



(17) 
(18) 



where M is as in (15) and (16). In addition we may assume that (14) holds with the term 
j||V 2 e||| 2 ( J 4 r ) instead of ||| V 2 e|| 2 2 t Ar j on the right hand side. Then the same argument as in 
(16) yields in addition that 



2^ 



\B 2 - M 



V 2 e\ z dx < din ( In - ) + C 4 



(19) 



These additional estimates can be used to improve the prcfactor of In (in j^j in the lower bound 
from 1 to i + e. 
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